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Measurement-based quantum computation (MQC) is a paradigm for studying quantum computation using many-body entan¬ 
glement and single-qubit measurements. While MQC has inspired wide-ranging discoveries throughout quantum information, 
our understanding of the general principles underlying MQC seems to be biased by its historical reliance upon the archetypal 
2D cluster state. Here, we utilize recent advances in the subject of symmetry-protected topological order (SPTO) to intro¬ 
duce a novel MQC resource state, whose physical and computational behavior differs fundamentally from the cluster state. 
We show that, in sharp contrast to the cluster state, our state enables universal quantum computation using only measure¬ 
ments of single-qubit Pauli X, Y, and Z operators. This novel computational feature is related to the “genuine” 2D SPTO 
possessed by our state, and which is absent in the cluster state. Our concrete connection between the latent computational 
complexity of many-body systems and macroscopic quantum orders may find applications in quantum many-body simulation 
for benchmarking classically intractable complexity. 


I. INTRODUCTION 

The idea of measurement-based quantum computa¬ 
tion (MQC), where computation is carried out solely 
through single-qubit measurements on a fixed many-body 
resource state and classical feed-forward of measurement 
outcomes m, is quite surprising. This is because it 
highlights the origins of quantum advantage in terms 
of entanglement and non-commutative measurements, 
uniquely quantum effects without counterparts in classi¬ 
cal mechanics. In particular, so-called universal resource 
states, the states that are capable of efficiently imple¬ 
menting universal MQC, represent a class of maximal 
entanglement in the classification of many-body entan¬ 
glement 0 , so that the structure and complexity of their 
entanglement is of great interest for advancing the under¬ 
standing of quantum computation. Following the canon¬ 
ical example of the 2D cluster state [5], many other uni¬ 
versal resource states have been found, including cluster 
states defined on various lattices [I], some tensor network 
states EBEflm, and model ground states in condensed 
matter physics such as 2D Affleck-Kennedy-Lieb-Tasaki 
(AKLT) states [TD1 - 1T5] . 

Given the existence of these various known universal 
resource states, a natural question is whether we might 
be able to find any common key feature, so as to explore 
more their variety in fundamental structures as well as 
practical applications. While the earliest resource states 
for MQC were found in short-range correlated states de¬ 
scribed as somewhat artificial tensor network states 013 - 
1T0] . a new insight has been that a class of short-ranged 
entangled states structured by symmetry, endowed with 
so-called symmetry-protected topological order (SPTO) 
riM24| , make excellent candidate resource states system¬ 
atically. Indeed, in the setting of ID spin chains, the 
ground states of several SPTO phases have already been 


shown to possess entanglement which can be leveraged 
to achieve various quantum computational tasks [3SH2I] ■ 
Here, in adopting the concept of SPTO, we carry out 
such an investigation for the first time in 2D MQC, and 
discover a completely new kind of MQC universal re¬ 
source state. Specifically, we first examine the 2D clus¬ 
ter state as well as a wide range of other universal re¬ 
source states, and show that their 2D SPTO is trivial, of 
the same nature as unentangled product states. Look¬ 
ing more closely, we find that these previously known 
resource states do possess some “weaker” SPTO, but es¬ 
sentially of a type closer to that of ID spin chains. Our 
discovery is made possible owing to the recent progress 
of research into SPTO, which has revealed a hierarchy 
of SPTO as representing different levels of nonlocality of 
quantum information (see the next section for details). 
We then introduce our “Union Jack” state, which in con¬ 
trast possesses SPTO entirely of a 2D nature, and demon¬ 
strate that it is not only a universal resource state but 
additionally is “Pauli universal,” in that it can imple¬ 
ment arbitrary quantum computation using only single¬ 
qubit measurements in the Pauli bases. As elaborated 
later, this feature is forbidden in the 2D cluster state 
on account of the Gottesman-Knill theorem [32], which 
proves the efficient classical simulability of certain quan¬ 
tum gates. We will conclude with the outlook that our 
proof of principle result about Pauli universality may be 
true for more general 2D SPTO resource states, which 
we connect to a possible deep connection between a hi¬ 
erarchy of SPTO in condensed matter physics and the 
so-called Clifford hierarchy of quantum computation. 


II. BACKGROUND 

A. Measurement-based Quantum Computation 
and the Clifford Hierarchy 
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Measurement-based quantum computation (MQC) is a 
means of utilizing an entangled resource state to perform 
computation using (generally adaptive) single-qubit mea- 




surements. Given a particular resource state, we specify 
our computational process by choosing a specific pattern 
of single-qubit measurements. Due to the probabilistic 
nature of measurement, different measurement outcomes 
will generally implement different computations. How¬ 
ever, rather than attempting to correct for unintended 
measurement outcomes at every step, we can instead 
represent the effect of such outcomes as the product of 
our intended operation with a so-called byproduct oper¬ 
ator. When these byproduct operators are sufficiently 
simple (e.g. Pauli operators), we can commute them 
through much of our computation, allowing disjoint mea¬ 
surements to be performed in parallel without adaptation 
of our measurement settings. 

The canonical MQC resource state is the 2D cluster 
state [3], which is a universal resource state, in that arbi¬ 
trary quantum circuits can be simulated efficiently using 
an appropriate sequence of arbitrary single-spin measure¬ 
ments m- The 2D cluster state is formed by preparing 
qubit states |+X) = ^(| 0 ) + | 1 )) on the vertices of a 
square lattice (with open boundary conditions), and ap¬ 
plying entangling controlled-Z ( CZ ) operations, defined 
in the computational basis by CZ \a,fi) = (—l) a/J | a, (3), 
between nearest-neighbor qubits. It is described by sta¬ 
bilizer generators, 

sW=xW 0 Z u \ ( 1 ) 

j'eneigh(i) 

where neigh(i) is the set of nearest neighbors of site i. An 
n-qubit cluster state | ipc) is the unique state satisfying 

s c ] IV’c) = IV’c) for i = 1,2,..., n. 

The Clifford hierarchy is an ordered collection of uni¬ 
tary gates of increasing computational generality [55] . 
The unitary gates in the d’th level of the Clifford hi¬ 
erarchy Cd are defined inductively, with C\ consisting of 
tensor products of Pauli operators, and Cd+i = {U\ VPs 
C±,UPW C Cd}- Each level of the Clifford hierarchy 
represents a greater degree of quantum-gate complexity 
in that, intuitively speaking, higher levels contain gates 
which are more “quantum” than those in lower levels. 
The gates in C 2 form a group, known as the Clifford 
group, which preserves the group of Pauli operators un¬ 
der conjugation. Exploiting this fact, the Gottesman- 
Knill theorem ;32] gives an efficient means of classically 
simulating any poly-sized circuit composed of gates in 
C 2 > provided that initialization and measurement occur 
in the single-qubit Pauli bases. By contrast, the gates in 
C 3 form a universal gate set for quantum computation. 

In MQC, a stronger notion of universality for resource 
states is Pauli universality, where the measurements used 
to carry out MQC are only of single-qubit Pauli opera¬ 
tors A', Y, or Z. While the 2D cluster state is a uni¬ 
versal resource state, it is formed from CZ gates in C 2 
and therefore can be efficiently classically simulated when 
only Pauli measurements are used, making the cluster 
state not Pauli universal. 


B. Symmetry-Protected Topological Order 

Symmetry-protected topological order (SPTO) fTfHT2H~) 
is a many-body phenomenon arising from many-body en¬ 
tanglement present in quantum states invariant under a 
symmetry group G. Given a state defined in d spatial 
dimensions with a finite correlation length, we say that 
this state has nontrivial d-dimensional SPTO precisely 
when it cannot be reduced to a product state using a 
finite-depth quantum circuit whose gates are of constant 
size and commute with G. In this sense, nontrivial SPTO 
can be thought of as an indicator of persistent entangle¬ 
ment, protected by G. More generally, two d-dimensional 
states are said to be in different (d-dimensional) SPTO 
phases when they cannot be transformed into each other 
using such a finite-depth, symmetry-respecting quantum 
circuit. 

Mathematically, d-dimensional SPTO phases are clas¬ 
sified by elements of T~L d+1 (G, 17(1)), the (d + l)’th coho¬ 
mology group of G, with the identity element of the group 
corresponding to the trivial phase of G-invariant product 
states (see Appendix [A] for an introduction to group co¬ 
homology theory). For example, when G = Z 2 there is 
only one (trivial) ID SPTO phase, but there are two 2D 
SPTO phases, one trivial and one nontrivial. Nontrivial 
SPTO can be detected and characterized by examining 
the manner in which G acts on edge degrees of freedom 
when a state is prepared on a manifold with boundaries 
[34H3Z]. Nontrivial ID SPTO manifests as a product 
of individual “fractionalized” degrees of freedom on the 
edge, which transform under projective representations 
of G. On the other hand, nontrivial 2D SPTO manifests 
in the form of long-range correlated edge modes, which 
transform under non-separable matrix product unitary 
representations of G [53] . Concrete examples of this dis¬ 
tinctive behavior of ID and 2D SPTO are shown in Fig¬ 
ure m 

An important—and often neglected—fact is that states 
in d spatial dimensions can be classified not only by a la¬ 
bel specifying its d-dimensional SPTO phase, but also 
by other labels associated with fc-dimensional SPTO, 
for 0 < k < d [T8] . We call this collection of 

SPTO labels the SPTO signature of a state, denoted 
by fid in d dimensions. For d = 2, H 2 has the form 
= ((© 2 ; B^,©^; ©o)), with ©& denoting a fc- 
dimensional SPTO label. For general d, fid contains 
( d ) fc-dimensional SPTO labels, corresponding to the 
(^) independent fc-dimensional surfaces in d-dimensional 
space. When classifying phases, the ©*, labels are cho¬ 
sen from 'H fc+ 1 (G, [/(l)), the collection of fc-dimensional 
SPTO phases for symmetry G. However, since we are 
concerned here mainly with the existence of nontriv¬ 
ial SPTO, we will use an abbreviated notation where 
Ok = 0 or 1 indicates trivial or nontrivial fc-dimensional 
SPTO, respectively. Unlike d-dimensional labels, the 
lower-dimensional components of a state’s SPTO signa¬ 
ture can be altered by a local G-symmetric quantum cir- 
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FIG. 1. Manifestation of ID and 2D SPTO in boundary sym¬ 
metry operators, where A', Y, Z, and CZ represent the appli¬ 
cation of the corresponding unitary operation. The transver¬ 
sal application of X is only a symmetry when each state is 
prepared on closed boundaries. Near edges of the system, 
the symmetry operator must be augmented with additional 
boundary terms, which reflect the distinct nature of ID vs. 
2D SPTO. a) The cluster state is invariant under the applica¬ 
tion of X to all sites within a region with closed boundaries. 
When X is instead applied to a region with open boundaries 
(boxed area), we must add additional Z (and hence Y. when¬ 
ever X and Z overlap) gates near the edges to obtain a genuine 
symmetry, b) The Union Jack state is invariant under the ap¬ 
plication of A to a region with closed boundaries. When X 
is instead applied to a region with open boundaries (boxed 
area), we must add additional CZ gates near the edges to 
obtain a genuine symmetry. The higher-dimensional SPTO 
manifests here as a symmetry representation which doesn’t 
factorize into disjoint unitaries, and is built from gates at a 
higher level of the Clifford hierarchy. 


cuit. However, these labels are unchanged by quantum 
circuits which respect both on-site and lattice transla¬ 
tional symmetries. See Appendix [A] for details about 
SPTO signatures. 


III. TRIVIAL 2D SPTO OF THE 2D CLUSTER 
STATE 


In this section, we determine the SPTO signature of 
the 2D cluster state, stated in Theorem [T] 

Theorem 1. The SPTO signature of the 2D cluster 
state with respect to on-site (Z 2) 4 symmetry is = 

((0; 1,1; 0)), corresponding to trivial 2D SPTO and 
nontrivial ID SPTO. 

The on-site (Z 2) 4 symmetry of the cluster state comes 
from treating a 2 x 2 unit cell as a single site, as shown in 
Figure [2ji. We refer to the four qubits within a unit cell 
by the labels NW, NE, SE, and SW. From Eq. |l]), we 
see that the global application of X to any of these four 
classes of qubits preserves the cluster state stabilizers, 
giving the system (Z 2) 4 on-site symmetry. This is the 
largest on-site symmetry group of the cluster state, and 
its SPTO phase with respect to this group sets its SPTO 
phase with respect to any on-site symmetry subgroup 

m- 
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FIG. 2. a) Part of the 2D cluster state on a square lattice, 
with 2x2 unit cells shown. The four generators of the (Z 2) 4 
on-site symmetry are labeled, b) Part of the circuit which 
disentangles the 2D cluster state. Solid lines indicate a CZ 
applied between two sites. The gate Ve is shown in center, 
which is the product of 6 CZ operations between sites (4, 8), 
(8,12), (12,14), (14,10), (10,6), and (6,4). Also shown are 
portions of the Ve gates directly above and below. Due to the 
“diagonal” CZ 's of adjacent Ve ’s canceling, a global tiling 
of these gates applies CZ between all adjacent NE and SE 
sites. This tiling is done in two layers, so that the gates in 
each layer don’t overlap. By applying displaced and rotated 
versions of these gates, we arrive at a symmetry-respecting 
circuit of depth 8, which disentangles the 2D cluster state to 
a trivial product state. 




We prove the 2D part of Theorem [T] by constructing a 
finite-depth quantum circuit, shown in Figure (2jo, whose 
gates each respect the on-site symmetry of the cluster 
state, but which disentangles the state to a trivial prod¬ 
uct state. Because the 2D component of a state’s SPTO 
signature is invariant under local symmetric quantum cir¬ 
cuits [181 , this suffices to prove our claim. A more careful 
analysis of the 2D cluster state is needed in order to prove 
its nontrivial ID SPTO. In Appendix [C] we study a pro¬ 
jected entangled pair state (PEPS) [333] representation of 
the cluster state, which lets us characterize the transfor¬ 
mation of its boundary under the (Z 2) 4 symmetry [37] . 
We find that individual sites along both horizontal and 
vertical boundaries transform under a projective repre¬ 
sentation of (Z 2 ) 4 , giving us a “smoking gun” indication 
of nontrivial ID SPTO. This fact, demonstrated rigor¬ 
ously in Appendix [C] completes our proof of Theorem [l] 


Importantly, a similar analysis of edge modes can be 
used to prove results analogous to Theorem [l] for many 
other known universal resource states, including clus¬ 
ter states defined on various lattices [3] and certain 2D 
AKLT states fT0Hl5| . In this sense, the impact of The¬ 
orem |T] is that not just the cluster state, but in fact the 
majority of commonly studied universal resource states, 
are characterized by the absence of 2D SPTO, with at 
most nontrivial ID SPTO. 
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FIG. 3. a) The Union Jack lattice on which our resource 
state is defined. Every vertex represents a qubit initialized in 
a |+A) state, and every triangular cell represents an applied 
3-body unitary CCZ. A 2x2 unit cell is shown, with respect 
to which our system has (Z 2) 3 symmetry generated by X 
applied to sites a, b, or c. The Z 2 symmetry of this state is 
a subgroup of (Z 2) 3 generated by applying A' to all sites, b) 
Measuring the control sites (red) in the computational basis 
collapses the remaining system into a random graph state. 
The edges of the graph lie on the “domain walls” between 
different control site outcomes. 


IV. THE RESOURCE STATE WITH NONTRIV¬ 
IAL 2D SPTO 

In this section we present a new MQC resource state 
that is both Pauli universal and possesses nontrivial 2D 
SPTO, as summarized in Theorem [2j This is in contrast 
to the 2D cluster state, which is universal but not Pauli 
universal, and only possesses ID SPTO. Our “Union 
Jack” resource state is composed of qubits, each of which 
is located at a vertex of the Union Jack lattice shown 
in Figure [3^,. It is constructed by preparing a |+X) 
state at every vertex, and then applying a 3-body doubly 
controlled-Z unitary operation, CCZ, to every triangu¬ 
lar cell in the lattice. CCZ is diagonal in the qubits’ 
computational basis with non-zero matrix elements: 


(* 1 * 2 * 3 ! CCZ \i x i 2 if) 


-1, if (* 1 , * 2 ,* 3 ) = (1,1,1), 
+1, otherwise, 


( 2 ) 

and belongs to the 3rd level of the Clifford hierarchy C 3 . 
The stabilizers generated by these gates are 


S^j=X^ 0 CZ^’ k \ (3) 

Cifc)etri(i) 

where (j, k ) £ tri(i) refers to all pairs of sites (j, k) which, 
together with i, form a triangle in the lattice of Figure[3[r. 
Our resource state | ifuj) is the unique state satisfying 
SujI^uj) = | 4>uj) for i = 1,2,...,n. Note, however, 
that it is not a so-called stabilizer state because its sta¬ 
bilizer group is not contained in the n-qubit Pauli group. 

Our resource state is an example of a “renormaliza¬ 
tion group (RG) fixed point” state used previously to 



Cluster 

Union Jack 


State 

State 

SPTO 

ID 

2D 

Formation 

c 2 

c 3 

Circuit 

Byproduct 

Operators 

Ci 

c 2 

Universal 

c 2 

Ci 

Measurements 


TABLE I. A summary of the SPTO present in our repre¬ 
sentative resource states, the quantum circuit used to form 
each state, the logical byproduct operators appearing during 
a computation, and the single-qubit operators whose eigen- 
basis we need to measure to achieve universal MQC. Cd refers 
to gates chosen from the d’th level of the Clifford hierar¬ 
chy. Higher-dimensional SPTO is associated with a higher 
gate complexity in the formation circuit and logical byprod¬ 
uct operators, and consequently requires less complexity to 
be added in the form of measurements. By contrast, we must 
perform measurements in eigenbases of gates from C 2 in order 
to achieve universal MQC with the cluster state. 


study properties of Z 2 SPTO [18j , and consequently has 
Z 2 symmetry. However, if we redefine a single site of 
our system to be a particular 2x2 unit cell (shown in 
Figure [3^,), then our system in fact has (Z 2 ) 3 symme¬ 
try. With respect to this latter group, our resource state 
can be seen as an example of a d = 2 decorated domain 
wall (DDW) state [30] , a method for creating systems 
with d-dimensional Z 2 x G SPTO in terms of systems 
with (d— l)-dimensional G SPTO (here G = (Z 2 ) 2 ). We 
should however emphasize the importance of our state 
being defined on the Union Jack lattice for proving The¬ 
orem^ as the 2D state in jT8j[40] is essentially defined on 
a triangular lattice, so that it disallows the intersection 
of domain walls under the procedure we use below for 
locally converting to a graph state mm, and thus may 
not be a universal resource state. On the other hand, 
our state is also an example of a generalization of graph 
states, called hypergraph states in the quantum informa¬ 
tion community mm, although their application for 
MQC has not been studied previously. 

Theorem 2. The Union Jack state is a Pauli universal 
resource state for MQC, meaning that arbitrary quantum 
circuits can be efficiently simulated using only measure¬ 
ments of single-qubit Pauli operators and feed-forward of 
measurement outcomes. Furthermore, its SPTO signa¬ 
ture with respect to on-site (Z 2 ) 3 symmetry is = 

(( 1; 0,0; 0)), corresponding to nontrivial 2D SPTO and 
trivial ID SPTO. 

Note that while we phrase Theorem [2] in terms of our 
state’s (Z 2 ) 3 SPTO, the same statement holds if we re¬ 
place (Z 2 ) 3 by Z 2 . Here we demonstrate the Pauli uni¬ 
versality of our state by efficiently simulating quantum 
circuits composed of Hadamard ( H ) and Toffoli (TOFF) 
gates — a universal set of unitary gates [13 — using only 
measurements of single-qubit Pauli operators. 
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FIG. 4. a) Our “interaction gadget”, which implements the non-Clifford operation Ui , and is formed by measuring X and 
Z on seven logical sites, Y on one control site, and Z on the surrounding control sites. Postselection on 13 of the latter Z 
measurement outcomes is required in order for us to connect our gadget to the surrounding cluster region, however this only 
introduces a constant overhead to the number of sites measured in our protocol, as shown in Appendix [D] b) A gadget for 
implementing SWAP within a cluster region. This allows us to implement nonplanar wire crossings, which are necessary for 
simulating arbitrary circuits composed of H and TOFF gates, c) A protocol for implementing a logical CCZ gate, where 
solid lines indicate teleportation of logical qubits. The majority of the sites involved have been converted to an extended 
cluster region, with the exception of the sites used to construct interaction gadgets. Our diagram only reflects the topology of 
the relevant logical connections, whereas a realistic implementation would involve a detailed measurement pattern to perform 
teleportation throughout the cluster region, as well as a significantly greater distance between neighboring gadgets. More 
explicit details of our protocol can be found in Appendix [D| 


Our means of simulating these circuits using the Union 
Jack state are divided into two parts. We first show 
that portions of our state can be converted into “clus¬ 
ter regions”, regions which are locally identical to the 2D 
cluster state. These cluster regions are used to prepare 
and readout qubit states, teleport these states over ar¬ 
bitrarily long distances, and apply Clifford gates (which 
include H gates) to them, all with the use of only Pauli 
measurements. We then demonstrate that we can imple¬ 
ment CCZ using certain “interaction gadgets”, which are 
prepared using Pauli measurements. Since we can imple¬ 
ment both H and CCZ gates, and because TOFF and 
CCZ are related by TOFF < 123 ) = H^CCZ^H^, 
the combination of cluster regions and interaction gad¬ 
gets lets us implement H and TOFF gates, and therefore 
arbitrary quantum circuits. 

Our technique for creating cluster regions within the 
Union Jack state is to induce a symmetry-breaking phase 
transition from 2D to ID SPTO. This involves first per¬ 
forming a computational basis measurement of all the 
Union Jack control qubits, shown in Figure [3jo. This 
symmetry-breaking measurement forces the remaining 
part of our system, which lives on a regular square lat¬ 
tice, into a random graph state whose edges (associated 
with nontrivial ID SPTO) appear along the domain walls 
in our measurement outcomes. In particular, we obtain 


an edge (CZ gate) in our graph whenever two adjacent 
measurement outcomes differ, and no edge whenever they 
agree. We can then use the exact same protocol as in 
[35] to reduce this random graph state to a state which 
is locally identical to the regular 2D cluster state. This 
protocol succeeds with a probability that converges expo¬ 
nentially fast to either 0 or 1 in the limit of large cluster 
regions, depending on whether our random graph state 
percolates and has a macroscopic spanning cluster of con¬ 
nected vertices. We perform numerical simulations of this 
percolation problem for different system sizes, and con¬ 
clude (see Figure [5]) that our Union Jack system is in a 
supercritical percolation phase and thus can be used to 
efficiently prepare connected cluster regions. 

Our technique for preparing interaction gadgets in¬ 
volves taking a small area of the Union Jack state 
and applying an appropriate pattern of Pauli mea¬ 
surements to it (see Figure Bo- When embedded 
within a cluster region, these gadgets implement a 
three-body non-Clifford logical gate Ui, defined by 

t® 23) = CCZ^ 12 ^ VCZ I ' 12) VCZ < ' 23) , where ^fCZ acts 
as \/CZ \a,P) = {i) a P \a, (3). Using [//, we can ob¬ 
tain CCZ by applying Ui three times to the same 
triple of qubits, but with the qubits cyclically per¬ 
muted each time. This permutation involves cross- 
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FIG. 5. A simulation of our percolation problem with increas¬ 
ing linear size, L. The exponential decay of the non-spanning 
probability is characteristic of the percolation supercritical 
phase, demonstrating that portions of our Union Jack state 
can be locally reduced to a 2D cluster state with arbitrarily 
high probability. These cluster regions are used to perform 
Clifford operations upon our computational qubits, as well as 
to shuttle these qubits between spatially separated interaction 
gadgets, which can be connected together to produce logical 
CCZ gates. 


ing adjacent wires, something which is forbidden in a 
strictly planar graph structure, but we can simulate 
a nonplanar wire crossing using a SWAP operation 
within our cluster regions (see Figure [4 Jd) . The iden¬ 
tity t/f 23) 17 f 31) [/f I 1 2 ) = CCZA^CZ^CZ^CZ^ 23 ') 
shows that this gives the desired operation of CCZ , up 
to byproduct CZ gates. These byproduct gates, as well 
as other byproduct Clifford gates which appear in our 
protocol, are adaptively eliminated within cluster regions 
by applying the appropriate inverse Clifford operations. 
This adaptive cancellation of byproduct operators is gen¬ 
erally necessary before the application of subsequent H 
or CCZ logical gates, since attempting to commute them 
through these gates would lead to a byproduct group 
which doesn’t close at any level of the Clifford hierarchy. 
Additional information about our protocol for establish¬ 
ing Pauli universality of the Union Jack state is given in 
Appendix [D] 

Looking at the proof just given, we see that the dis¬ 
parity between the universality of the cluster state and 
the Pauli universality of the Union Jack state arises from 
the difference in gates implementable by each state un¬ 
der Pauli measurements, C 2 for the former and C 3 for 
the latter. Some insight can be gained by comparing 
this computational difference to the fact that the clus¬ 
ter state and Union Jack state respectively possess ID 
and 2D SPTO, as summarized in Table [T] Generalizing 
from these examples, we might expect this correspon¬ 


dence between SPTO and the Clifford hierarchy to ex¬ 
tend to a wider class of SPTO states, providing a general 
link between types of SPTO and degrees of quantum gate 
complexity. Such a correspondence was demonstrated in 
03 EH] for topological quantum error-correcting codes, 
but proving this in the setting of MQC would give a 
means of directly associating computational character¬ 
istics to SPTO states, without the need for an auxiliary 
higher-dimensional topologically ordered system. 

V. DISCUSSION 

Although pertaining most immediately to MQC, our 
main results can be fruitfully interpreted as general state¬ 
ments about the interplay of two intrinsically quantum 
ingredients, entanglement and measurement, which play 
a leading role in quantum information science. Our The¬ 
orem [T] demonstrates that previously studied resource 
states, despite differing in their microscopic details, pos¬ 
sess identical forms of macroscopic entanglement, namely 
ID SPTO. While such entanglement is sufficient for uni¬ 
versal quantum computation using arbitrary single-qubit 
measurements, our Theorem [2] demonstrates that the 
use of more complex forms of entanglement, namely 2D 
SPTO, lets us achieve the same results using simpler 
Pauli measurements. As argued in the previous section, 
we expect that this tradeoff between entanglement and 
measurement is not only true of more general quantum 
systems, but in fact is evidence of a deep connection 
between the hierarchies of SPTO and the Clifford hier¬ 
archy of quantum computation. Such a connection be¬ 
tween the computational complexity of many-body sys¬ 
tems and their emergent macroscopic behavior would 
give a means of converting canonical condensed matter 
tools, such as order parameters, into interesting indica¬ 
tors of computational behavior, as was done with ID 
spin chains in |30j . The natural connection we demon¬ 
strate between the computational complexity of many- 
body systems and emergent macroscopic order may find 
applications for better understanding the emergence of 
classically intractable complexity within quantum many- 

body simulation 03150]. 
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In Appendix [A] we provide a more complete discus¬ 
sion of SPTO, including an overview of the relevant con¬ 
cepts from group cohomology theory. In Appendix |B] we 
demonstrate the precise manner in which the 2D cluster 
state and the Union Jack state are examples of SPTO 
fixed-point states originally introduced in [15]. Finally, 
in Appendices [C] a nd [D] we give the full proofs of our 
Theorems [T] and|2| 


Appendix A: Symmetry-Protected Topological Or¬ 
der 

We give here a more complete discussion of SPTO, 
and in particular the possible SPTO signatures that are 
allowed for an arbitrary 2D state. We restrict our dis¬ 
cussion to systems with an on-site symmetry G, and ig¬ 
nore SPTO arising from global symmetries, such as time 
reversal, spatial inversion, or lattice point group symme¬ 
tries. However, we do consider the effect of lattice trans¬ 
lational symmetries, since this symmetry is necessary for 
lower-dimensional portions of our SPTO signature to be 
well-defined. After having given this general discussion of 
SPTO, we state the classification of several SPTO phases 
in 2D and ID which are relevant for our purposes. 

The classification of SPTO phases is closely tied to 
group cohomology theory, so we first give a brief in¬ 
troduction to some of the concepts from that field. 
Given a symmetry group G, we can construct n-cochains 
w n , which are functions from the direct product of n 
copies of G to the group of complex phases, [7(1) = 


{a £ C | aa* = 1}. The collection of n-cochains 
form an abelian group C n (G , [7(1)) under pointwise mul¬ 
tiplication, with the product of cochains co n and u>' n 
given by a cochain uj n u' n , where (u n w^)(gi,..., g n ) = 
U] n {gi, ■ • •, gnWnigi, ■ ■ ■, g n )- The identity element in 
C n (G , [7(1)) is the trivial n-cochain, co°(gi,... ,g n ) = 1. 
We define an operation called the coboundary operator, 
d„:C”(G,[7(l))^C" +1 (G,[7(l)), by 


(d n uj n )(gi ,..., g n + i) — 

w n (52, • • •, g n + i)o4 _1) " +1 (51,, g n ) 

n 

n ^ (flu • • • ,9k-l,9k9k+l,9k+2, ■ ■ ■ , g n +l)- (Al) 

k— 1 

A special role is played by the n-cocycles and n- 
coboundaries, which form subgroups of C n (G,U( 1)) de¬ 
noted by Z n (G,U( 1)) and B n (G,U( 1)), respectively. An 
n-cochain is an n-cocycle (resp. n-coboundary) if it 
lies in the kernel of d n (resp. the image of d n -\). 
More explicitly, Z n (G , [7(1)) = {u n \ d n uj n = w° +1 } and 
B n (G,U{ 1)) = {w„ | 3w n _i s.t. d n oj n -i = l o n }. One 
can show that the composite of coboundary operators 
d n and d n +1 is trivial, in that it sends every n-cochain 
to the identity (n + 2)-cochain. This implies that ev¬ 
ery n-coboundary is an n-cocycle, so that B n (G, [7(1)) C 
Z n (G,U( 1)). 

We define the n’th cohomology group of G, 
r H n (G,U( 1)), to be the (abelian group) quo¬ 
tient of Z n (G, [7(1)) with respect to B n (G,U( 1)), 
H n (G,U{ 1)) = Z n (G,U(l))/B n (G,U(l)). Equivalently, 
this is the group of equivalence classes of n-cocycles, 
TL n {G, [7(1)) = {[co n ]\uj n £ Z n (G,U(l))}, under the 
equivalence relation [uj n \ = [w(J O u„ = where 

w" is an arbitrary n-coboundary. For w n £ Z n (G,U( 1)), 
we will call [u> n ] £ r H n (G,U( 1)) the cohomology class 
associated to cu n . 

The relevance of this discussion for our purposes is that 
SPTO phases of G-invariant many-body systems living 
in d-dimensional space are classified by elements of the 
(d + l)’th cohomology group. In particular, it was shown 
in [T5] that given any two distinct cohomology classes in 
H (d+1 )(G,U{ 1)), we can construct d-dimensional “fixed 
point” systems labeled by the cohomology classes which 
belong to different SPTO phases. This construction is 
discussed in more detail in Section El 

An important point is that systems with both on¬ 
site G symmetry and translational symmetry admit a 
richer classification of SPTO phases |T5]. In particular, 
while the SPTO phase of a system without translational 
symmetry can be uniquely classified by a single coho¬ 
mology class, with additional translational symmetry in 
place, the SPTO phase is classified by a full SPTO sig¬ 
nature fid, which consists of an ordered list of different 
cohomology classes. For systems in 2D, this signature 
is of the form fl 2 = (([w 3 ]; [w^], [c*4^]; [aq])), with 
M G n 3 (G,U( 1)), [4 x) ],[4 y) ] G H 2 (G, [/(!)), and 


[wi] G H 1 {G,U(\)). We refer to these respectively as 
the 2D, ID, and OD portions of fl 2 . For SPTO sys¬ 
tems in d physical dimensions, there will generally be 
(^) components to the Axdimensional sector of the SPTO 
signature, corresponding to the number of independent k- 
dimensional surfaces in d-dimensional space. Due to our 
present focus on only whether or not a system possesses 
SPTO, we often use an abbreviated means of writing the 
components of an SPTO signature, wherein a phase label 
is written as 0 if it corresponds to the trivial phase, and 
as 1 if it corresponds to any nontrivial phase. 

We now introduce a few examples of concrete SPTO 
phases in 2D and ID associated with various symmetry 
groups. Since there is always a trivial phase for every 
symmetry group and dimension, we will often neglect to 
mention these phases. 

For G = Z 2 , we have no nontrivial phases in ID, and 
one nontrivial phase in 2D. Our Union Jack state lives in 
this nontrivial 2D Z 2 phase when its symmetry group is 
taken to be Z 2 . 

For G = D 2 — (Z 2 ) 2 , we have one nontrivial phase in 
ID (known as the D 2 Haldane phase), and 7 nontrivial 
phases in 2D. D 2 is the smallest symmetry group which 
is capable of manifesting SPTO in ID. 

For G = (Z 2 ) 3 , we have 7 nontrivial phases in ID 
and 127 nontrivial phases in 2D. Using a known de¬ 
composition of 2D abelian SPTO phases (those with 
G abelian), we can structure the 2D (Z 2 ) 3 phases as 
H 3 ((Z 2 ) 3 ,[/(1)) ~ (Z 2 ) 3 x (Z 2 ) 3 x Z 2 [EE|. The first 
(resp. second) (Z 2 ) 3 factor encodes the “type I” (resp. 
“type II”) phases, those whose nontrivial SPTO arises 
from only one (resp., from pairs) of the Z 2 compo¬ 
nents in (Z 2 ) 3 . The last Z 2 in the decomposition of 
■H 3 ((Z 2 ) 3 , U(l)) is the unique “type III” component of 
the phase, which is due to a nontrivial combination of all 
three Z 2 components in (Z 2 ) 3 . Our Union Jack state with 
(Z 2 ) 3 symmetry belongs to the phase (0,0,1), meaning 
the unique phase with trivial type I and II SPTO, and 
nontrivial type III SPTO. 


Appendix B: The Union Jack and Cluster States as 
SPTO Fixed Point States 

In this Section, we demonstrate how both the Union 
Jack and 2D cluster states are examples of the construc¬ 
tion of H 8 for constructing special RG fixed point states 
with nontrivial SPTO from nontrivial cocycles of a sym¬ 
metry group G. We show how our Union Jack state be¬ 
longs to this class of states both for G = Z 2 and for 
G = (Z 2 ) 3 , and how the 2D cluster state belongs to this 
class of states. 

The construction of TS] gives a means of taking d- 
dimensional SPTO signatures, along with a representa¬ 
tive (k + l)-cocycle for each fc-dimensional component 
of the signature, and constructing a d-dimensional state 
with that SPTO signature. For our purposes, we will fo¬ 
cus on d = 2, for which the 2D, ID, and/or OD labels are 


allowed to be nontrivial. We will restrict first to the case 
of trivial lower-dimensional SPTO (the case considered 
almost exclusively in T8]), and later explain how these 
lower-dimensional labels can be made nontrivial. 

To construct a 2D state from a chosen group G and 
3-cocycle W 3 , we first choose a triangulated 2D lattice on 
which our state will live, and assign a Hilbert space Hq 
to every lattice vertex. Hq has dimension |G|, the order 
of G, and is spanned by an orthonormal basis labeled 
by the elements of G, {|.< 7 )} s eG- G acts on Hq as the 
regular representation uq , with u g \h) = \gh) for every 
g,h G G. We first initialize every Hq in the unique in¬ 
variant state |</>g) = (1/y/fGf) J2 g eG Id)’ which gives a 
symmetric global product state with trivial SPTO. We 
then apply to this system a collection of 3-body unitary 
gates, each formed from our chosen 3-cocycle, which gen¬ 
erates the nontrivial 2D SPTO. The 3-body unitary 013 
generated from a 3-cocycle W 3 is diagonal in the G-basis, 
and has non-zero matrix elements of 


(ghf\co 3 \ghf) =u: 3 (g,g 1 h,h 1 f). 


(Bl) 


Our desired state is obtained by applying u >3 or its in¬ 
verse to the vertices around every triangular cell in our 
chosen lattice. Whether we apply C 03 or W 3 to a particular 
triangular cell, as well as how we match up the 3 indices 
in Eq. (Bl) with the three sites around that cell, depend 


on a certain ordering of lattice vertices. While the full 
details are given in |TS] , if we restrict to 3-colorable lat¬ 
tices we can always choose each of the three indices to 
be matched up with a different vertex color in a fixed 
manner. 

Choosing G = Z 2 ~ {0,1}, this construction outputs 
qubit states, with \4>g) = |+AT). To produce our Union 
Jack state, we work with the Union Jack lattice, and 
choose our 3-cocycle to be 


( h j -1 ’ if (9,Kf) = (i,M) 

otherwise. <B2) 


Although this 3-cocycle produces a unitary W 3 which 
is distinct from CCZ, the global state it produces is 
nonetheless the same. This can be seen from the rela¬ 
tion d 4 123) = CCZ^CZ^\ which allows us to show 
that the transversal application of 0)3 to qubits in any 
3-colorable lattice with closed (nonexistent) boundary 
yields the same global unitary as the transversal applica¬ 
tion of CCZ. This proves that the Union Jack state is a 
Z 2 SPTO fixed point state, associated with the cocycle 
of Eq. (B2). Because this 3-cocycle belongs to the unique 
nontrivial cohomology class in "H 3 (Z 2 , t/(l)), our Union 
Jack state consequently has nontrivial 2D SPTO. 

Showing that our Union Jack state is isomorphic to a 
(Z 2 ) 3 SPTO fixed point state is less obvious, since the lat¬ 
tice vertices of such states aren’t associated with qubits, 
but rather with 8 -dimensional qudits. We can get around 
this difficulty by first treating each of the Z 2 factors in 
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(Z 2 ) 3 as a separate qubit system, and imagining these 
three factors to be stacked vertically in three layers at 
each lattice site. Note that this stacking is merely a 
convenient means of visualizing the separate qubit fac¬ 
tors in (Z 2 ) 3 , while our lattice remains a genuine 2D lat¬ 
tice. In this case, the state we initialize each site in is 
| <j>c) = |+A)® 3 , a tensor product of one |+A) state on 
each layer. If we write a generic element g £ (Z 2 ) 3 as 
g = (5ii 52 , 53 ), where each £ Z 2 is associated with the 
V th layer, then we can choose the following 3-cocycle 


^3 (g,h,f) = 


- 1 , if (gi,h 2 ,f 3 ) = ( 1 , 1 , 1 ) 
+ 1 , otherwise, 


(B3) 


where addition is modulo 2. Using the relation Eq. (B1 1 , 


we can show that ui 3 equals a CCZ gate on the qubits 
indexed by gi, /i 2 , and f 3 , along with other terms which 
cancel when uj 3 is applied globally. In other words, Cj 3 
ends up having a nontrivial action only on the qubits 
on the first layer of the first site acted on, the second 
layer of the second site, and the third layer of the third 
site. If we apply ui' 3 transversally to all triangular cells 
on a 3-colorable lattice, then at each site only one of the 
three layers is acted on nontrivially, with the other two 
layers remaining unchanged. Thus, using i2i' 3 to construct 
a (Z 2 ) 3 SPTO fixed point state defined on a Union Jack 
lattice with n vertices yields a state which is a tensor 
product of our Union Jack state on n qubits, with |+X) 
on the remaining 2 n qubits. This proves that, up to 
addition/removal of ancilla |+X) states, the Union Jack 
state is a (Z 2 ) 3 SPT O fi xed point state, associated with 
the cocycle of Eq. (B3l. This cocycle belongs to the 
nontrivial (Z 2 ) 3 cohomology class described at the end 
of Section [Aj which consequently specifies the nontrivial 
(Z 2 ) 3 SPTO phase our Union Jack state belongs to. 

As the 2D cluster state only possesses lower¬ 
dimensional SPTO, we must use an extended version of 
the previous construction to obtain the cluster state as 
an SPTO fixed point state. In m it is shown that to 
generate 2D fixed point states with ID SPTO, we can 
use a construction almost identical to that given above, 
but instead of starting with a 3-cocycle ui 3 and convert¬ 
ing it into a 3-body gate Cj 3 , we start with a 2-cocycle w 2 
and convert it into a 2 -body gate w 2 , which has non-zero 
matrix elements of 


{gh\ u 2 \gh) = u 2 (g,g 1 h). (B4) 

w 2 is then applied to all edges of a chosen 2D lattice, on 
which one copy of \4>g) has been prepared at every vertex. 
To generate the 2D cluster state in this manner, we can 
choose G = (Z 2 ) 2 and use a similar decomposition of the 
local Hilbert space into two qubits, stacked vertically in 
two layers. We then utilize the 2-cocycle 


W 2 (g,h) 


- 1 , if (gi,h 2 ) = ( 1 , 1 ) 
+ 1 , otherwise, 


(B5) 


where £ Z 2 is associated with the *’th component 

of g, h £ (Z 2 ) 2 . This 2-cocycle produces a 2-body unitary 
w 2 which upon global application is equivalent to a CZ 
gate on the qubits indexed by gi and /i 2 , and an identity 
gate on the rest of the qubits. In close analogy to how 
the Union Jack state was shown above to be a (Z 2 ) 3 
SPTO fixed point state, we can work with the 2-colorable 
square lattice and show that the transversal application 
of d > 2 to all edges of the lattice yields a state which is a 
tensor product of the 2D cluster state on n qubits, with 
|+A') on the remaining n qubits. This proves that, up to 
addition/removal of ancilla |+X) states, the cluster state 
is a (Z 2 ) 2 SPTO fixed point state. 

Finally, we note that some care is required regarding 
the symmetry group of the 2D cluster state. The con¬ 
struction we just outlined outputs the cluster state as 
an SPTO fixed point state with (Z 2 ) 2 symmetry, simi¬ 
lar to how the ID cluster state is most naturally seen as 
possessing nontrivial SPTO associated with (Z 2 ) 2 sym¬ 
metry. However, as seen from Eq. (C31, if we choose any 
particular (Z 2 ) 2 subgroup of the full (Z 2 ) 4 on-site sym¬ 
metry, we obtain a virtual representation of our symme¬ 
try which is non-projective in at least one direction. This 
leads to an SPTO signature which is either ((0 ; 0,1; 0)) 
or ((0; 1,0; 0)), rather than the SPTO signature of 
((0; 1,1; 0)) which appears in Theorem [lj We inter¬ 
pret this fact as an indicator that for states with lower- 
dimensional SPTO, we must take care in choosing the 
symmetry group we use to arrive at an SPTO signature. 


Appendix C: SPTO Signature of the 2D Cluster 
State 

We present here a full demonstration that the SPTO 
signature of the 2D cluster state is = ((0; 1,1; 0)), 
as stated in Theorem |T] To do this, we need to determine 
the various cohomology classes corresponding to different 
components of the cluster state’s signature. One known 
way 133 of doing this is by working with a projected 
entangled pair state (PEPS) description of the cluster 
state, and examining the behavior of the representation 
of its on-site symmetry group (Z 2 ) 4 along the boundary. 

Restricting to states which live on a square lattice, a 
PEPS representation consists of a rank-5 tensor, A £ 
H p (g> (H/)® 4 , where H p and H v are referred to as the 
physical and virtual Hilbert spaces, and where H* de¬ 
notes the Hilbert space dual to H. A can also be inter¬ 
preted as a map A : H* —► (H/)® 4 . We associate one 
copy of A to each site of our lattice, with H p correspond¬ 
ing to the Hilbert space of that site, and the four H*' s 
being used to represent correlations between our site and 
each of the four nearest-neighbor sites. The dimension of 
H v , D v , is the bond dimension of our PEPS representa¬ 
tion, and can be thought of as a measure of entanglement 
in the system. The condition for A to be a PEPS rep¬ 
resentation of a many-body state | ip) is that the “tensor 
trace” of the >4’s at every site, formed by contracting 
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a) 



FIG. 6. a) A single PEPS tensor for a square lattice. The dot¬ 
ted line represents our physical system, which corresponds to 
a single site of our lattice, and the four solid edges represent 
the virtual space, b) After assigning a PEPS tensor to every 
site of our lattice, we obtain a physical state by taking the 
“tensor trace” of all tensors. This involves contracting every 
pair of adjacent virtual indices using a maximally entangled 
state |(/>o) = Ef=”i M, with D v the virtual space dimen¬ 
sion. On a lattice with no boundary, this will contract out 
all of the virtual spaces, leaving only our physical many-body 
state c) An example of the physical/virtual symmetry 
correspondence given in Eq. ( |C1| ) for the 2D cluster state. 
Our PEPS tensor is defined relative to a 2 x 2 physical unit 
cell, with a four-qubit physical space and two-qubit virtual 
spaces. Different generators of (Z 2) 4 will produce different 
combinations of X and Z on the virtual space, whose non¬ 
commutativity demonstrates the nontrivial ID SPTO of the 
2D cluster state. 


every pair of adjacent H*’ s using maximally entangled 
states \cj)o) = |i, i), yields \ip). This condition is 

depicted in Figure [GJd. 

Given a PEPS representation A of our many-body 
state \ijj), the condition for \if>) to be invariant under 
our on-site symmetry G, whose physical representation 
is uq = { u g | g £ G}, is that there exists a virtual repre¬ 
sentation of G, Uq , such that 


Au G = e l9G U G A. (Cl) 

In other words, when A is seen as a map from the phys¬ 
ical to the virtual space, A is required to be (possibly up 
to phase) an intertwiner between the representations ug 
and Ug- e l9a = { e l9g \ g S G} is a unitary character of 
G, and using the fact that the collection of these charac¬ 
ters is isomorphic to'H 1 (G,/7(l)), the particular choice of 
e l9a ends up specifying the OD component of our SPTO 
signature. 

With the virtual representation U G '■ (H* )® 4 —► 
(ip*)® 4 in hand, we can calculate the remaining portions 


of the SPTO signature of our state \ip). The 2D portion 
of this signature relates to whether or not we can decom¬ 
pose Uq into a tensor product of four unitaries on the 
four virtual subsystems in (U*)® 4 . If we cannot, such 
that Ug is necessarily an entangled representation, then 
our state | ip) has nontrivial 2D SPTO. In such cases, 
there are several (somewhat involved) procedures for ex¬ 
tracting a 3-cohomology class to classify the 2D SPTO 
phase, but since our current interest is in the case of triv¬ 
ial 2D SPTO, we won’t discuss these here. The interested 
reader can consult El ED El] for examples of methods 
for obtaining information about 2D SPTO. 

Given trivial 2D SPTO, we can write Uq as a tensor 
product of four terms, which we will assume has the form 
Ug = 0 (U^)* 0 U^ ] 0 (U^y. These four terms 

correspond to, in order, the left, right, top, and bot¬ 
tom portions of our virtual representation, where ( Uq ■*)* 
(resp. (Uq' 1 )*) represent the complex-conjugated ver¬ 
sions of Uq ' 1 (resp. Uq' 1 ). We refer to Uq ' 1 and Uq^ 
as the horizontal and vertical components of our virtual 
representation, and these determine the ID portion of 
our SPTO signature. In particular, whether or not our 
system has nontrivial ID SPTO is equivalent to whether 
or not the horizontal/vertical components of our repre¬ 
sentation are nontrivial projective representations of G. 
More concretely, the product of two elements of U^\ 
Ug^ and ujy'* (g standing for either x or y), will gen¬ 
erally only equal U^ up to a phase factor, such that 
Ug^ujy ' 1 = <*4 fJ '\g,h)U^. Multiplication of elements 
of Uq 1> is associative, and this condition ends up forc¬ 
ing our phases uj^ig, h ) to be 2-cocycles. The coho¬ 
mology classes of these horizontal and vertical cocycles, 
[o;^] and then form the ID components of fl 2 , 

the SPTO signature of \ip). 

Let’s use these techniques to determine the SPTO sig¬ 
nature of the 2D cluster state. We can choose a PEPS 
representation for a single qubit site of the 2D cluster 
state as A I*) ® with the Ai € (P/**)® 4 
given by 


A 0 = <+X, 0, +X, 0| , A 0 = (-X, 1, -X,l\- (C2) 

H v is here a qubit space, and the ordering of our systems 
inEq. (C2) is as (H^ ,eft) ^H^ right) 0E< M 0^ 6otiom) )*. 
We are interested in the SPTO signature of the 2D clus¬ 
ter state with respect to a 2x2 unit cell, since the cluster 
state then has its maximal on-site symmetry group of 
G = (Z 2 ) 4 . To determine this, we c ontr act together four 
copies of the PEPS tensor of Eq. (C2) to form a 2 x 2 

PEPS tensor, Aq X2 \ and then find the virtual symme- 

(x) (v) 

try representations Uq and Uq. These each act on a 
two-qubit virtual space, which for Uq is decomposed 
as (Hv° pS> 0 £[( bottorn ')y ) anc j f or jj(v) j s decomposed as 
(Hi left) 0 Hi right) y. 
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As in the main text, we label the generators of (Z 2) 4 
by their respective locations in the 2x2 unit cell. One 
can then verify that the following choice of virtual sym- 

(2x2) 

metry representation makes our PEPS tensor A y c ' an 
intertwiner with respect to the physical representation 
uq (see Figure [6];): 


U { * ] w = Z®I 

= x ® I 

UgE = I® X 
Usw = I ® Z 


u { ^ w = Z®I 
Une = I ® Z 
= I® X 
uj& = X ® I 


(C3) 


The fact that we can choose a for m fo r Uq which fac¬ 
torizes into parts and satisfies Eq. (Cl) with e lSc = 1 is 
confirmation of the trivial 2D and OD SPTO of the 2D 
cluster state. The only thing that remains is determin¬ 
ing the two ID components of the SPTO signature. We 
can show that these are both nontrivial by considering 
the commutation relation of elements of Uq and Uq. 
Whi le ( Z 2) 4 is abelian, the virtual representations in 
Eq. (C3| aren’t, as shown by U^U^(U^) t (U^) t = 
U nwWw( U nwV( u swV = -I® 2 - This means that the 
2-cocycle associated with each of our virtual rep¬ 
resentations is different from the identity. Furthermore, 
multiplying either of these 2-cocycles by an arbitrary 2- 
coboundary is equivalent to modifying the phases asso¬ 
ciated to our individual Ug^ as Ug^ i-T uj\(g)Ug^\ with 
w\{g) € C 1 (G,U( 1)). This has no effect on the com¬ 
mutators of our symmetry group, which proves that our 

fx) (y) 

2-cocycles u > 2 and uv; are in nontrivial 2-cohomology 
classes. The SPTO signature of the 2D cluster state 
is therefore = ((0; 1,1; 0}), meaning trivial 2D 
SPTO and nontrivial ID SPTO, with the latter belong¬ 
ing to the nontrivial D 2 Haldane phase. 


Appendix D: Proof of the Pauli Universality of Our 
Resource State 

In this Section, we give a proof of the fact that our 
Union Jack resource state is Pauli universal, meaning 
that it can carry out universal MQC using only mea¬ 
surements of single-qubit Pauli operators. Achieving this 
universality requires several components, namely: 

• We can convert regions of our Union Jack to “clus¬ 
ter regions”, which are locally isomorphic to the 2D 
cluster state. This involves carrying out a pattern 
of computational basis measurements which con¬ 
verts (a part of) our state to a random graph state. 
The protocol of [3B] (which uses only Pauli mea¬ 
surements) is then used to concentrate this state 
into a 2D cluster state, which in turn requires the 
percolation problem associated with our random 
graph states to lie in a supercritical phase. We 



FIG. 7. A layout of our two-parameter percolation model. 
Cells labeled with pi (i = 1,2) are independently sampled, 
such that the probability of obtaining an outcome of 1 in 
that cell is pi. An edge of our random graph state is set 
when two adjacent nodes differ in their values. This yields 
a deterministically empty lattice at (pi,p 2 ) = ( 0 , 0 ) or ( 1 , 1 ), 
and a deterministically full lattice at (pi,P 2 ) = ( 0 , 1 ) or ( 1 , 0 ). 
Additionally, setting p\ = 0 (resp. P 2 = 0) gives a percolation 
problem which is isomorphic to a site percolation problem 
on a square lattice with a bond probability of P 2 (pi). Our 
problem of interest is located at (pi,P 2 ) = (§, 5 ). 


demonstrate the supercriticality of this percolation 
problem, and thereby the ability to prepare cluster 
regions within our state, in Section |D 1[ 

• We can teleport states and implement Clifford op¬ 
erations on them within the cluster regions of our 
state, using only Pauli measurements. Due to these 
cluster regions being identical to connected regions 
of the cluster state, we can use the same measure¬ 
ment patterns described in [2] to implement these 
Clifford operations, which use only Pauli measure¬ 
ments. 

• We can create “interaction gadgets”, which imple¬ 
ment a three-qubit non-Clifford operation, C/j 123 '* = 

CCZ^VCZ^ VCZ {23) , using only Pauli-basis 
measurements. Furthermore, these gadgets can be 
connected to a surrounding cluster region with a 
finite success probability, allowing us to use these 
gadgets as logical gates which we can connect to¬ 
gether to create a CCZ operation. We demonstrate 
these various facts in Section ID 21 

Taken together, these various facts successfully demon¬ 
strate the Pauli universality of our Union Jack state. 


1. Conversion to a 2D Cluster State 

After giving a more complete description of the reduc¬ 
tion of our Z 2 resource state to a random graph state, 
we describe the simulations we use to verify that the as¬ 
sociated percolation problem is indeed in the supercrit¬ 
ical phase. These simulations involve the construction 
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FIG. 8. The percolation phase diagram of our two-parameter 
model. Red (bottom left and upper right) indicates a subcrit- 
ical phase, while green (upper left to bottom right) indicates 
a supercritical phase. The yellow region contains the criti¬ 
cal line separating the phases. This division is based on the 
spanning probability p ap an when m = 100, and in particular 
whether p apan < 0.05, p apan > 0.95, or 0.05 < p span < 0.95. 
From the placement of our problem of interest at (pi,P 2 ) = 
(1/2,1/2), it is clear that we are within a supercritical phase, 
and can therefore use our 2D SPTO state as a universal re¬ 
source for MQC. 


of a two-parameter model which includes as a special 
case the percolation problem associated to our random 
graph state reduction protocol. We show that our partic¬ 
ular percolation problem lies within a supercritical phase, 
thus demonstrating that the protocol of [3B] can be used 
to efficiently convert these random graph states to a 2D 
cluster state with arbitrarily high probability. 

As described in the Methods, the method we use for 
reducing our 2D SPTO resource state to a random graph 
state consists simply of measuring all of the control sites 
in the computational basis. Given n control sites initially, 
upon measurement we obtain a string of random out¬ 
comes c = (ci, C 2 ,..., c n ). What is the reduced state of 
the logical portion of our system given a particular string 
of outcomes c? To figure this out, we exploit the fact that 
the projector associated with our measurement outcome 
commutes with all of the CCZ' s, since the latter are di¬ 
agonal in the computational basis. Thus, the state of our 
system after measurement is the same as if we had initial¬ 
ized the control sites in their post-measurement states, 
and afterwards applied CCZ everywhere in our lattice. 
The resulting (unnormalized) state is then 

tHc)) = -L n (cztyw+c'w \+x) 9n . (Di) 

e&c 2 

Here, C 2 is the collection of edges in our lattice, CZ t is 
a controlled-Z gate applied to the endpoints of a logical 


FIG. 9. a) The percolation probability for lattices of increas¬ 
ing linear size L, as we vary a parameter e from 0 to 1. The 
marginal bond probability varies as pb = e(l — |e), and the 
critical bond probability is seen to be ps = 0.484 ± 0.001. b) 
Using the same tools as were used in (a) to study the canoni¬ 
cal square lattice bond percolation problem. The critical bond 
probability is known to be 0 and our simulation reproduces 
this, locating it at pb = 0.500 ± 0.001. 


edge £, while c(£) and c'{t) are the measurement out¬ 
comes obtained on the two control sites adjacent to £. 
The factor of l/-\/2™ emerges from the inner product of 
our n measurement outcomes (0| or (1| with the |+X)’s 
which were used to initialize our state. What Eq. (Dl) 
tells us (ignoring normalization) is that whenever the 
measurement outcomes on two adjacent control sites are 
not equal, a CZ operation is performed on the logical 
edge in between them, while nothing is done when the 
measurement outcomes are the same. 

From this description, it is easy to see that every state 
■0(c)) is a graph state, whose edges lie only along do¬ 
main walls of the control site measurement outcomes. 
The control site outcomes themselves are uncorrelated 
and uniformly distributed, which follows from the equal 
mag nitu de of all of the unnormalized reduced states in 
Eq. (Dl). More precisely, the probability of obtaining a 
particular outcome c, p{ c), is given by 


p( C ) = (0(c) 0(C)) = —. 


(D2) 


Ignoring the quantum origin of the probabilities, this 
probabilistic reduction to a graph state can be seen as 
defining a (classical) percolation problem, wherein edges 
of a graph are filled based on the configuration of random 
control site variables. We wish to conclusively determine 
whether this percolation problem, with site probabilities 
given by Eq. (D2), corresponds to subcritical or super¬ 


critical behavior in the large-system limit. More explic¬ 
itly, from the known behavior of percolation problems, 
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we expect that the probability of obtaining a connected 
graph component which connects arbitrarily distant por¬ 
tions of our lattice goes to either 0 or 1 as we make our 
system size larger, and we would like to know which of 
these possibilities holds. 

To do this, we carry out numerical simulations of a 
two-parameter percolation model identical to ours, but 
with tunable probabilities for different control site out¬ 
comes. While Eq. (D2) corresponds to a probability of 
| of obtaining 1 on any arbitrary control site, our vari¬ 
able model has probabilities of p\ on one half of the sites, 
and P 2 on the other half of the sites. Figure [7] shows the 
checkerboard-style layout of these sites. The percolation 
problem defined by our actual system then corresponds 
to the point Pi = P 2 = 1/2. 

Figure [8] shows a phase diagram of this two-parameter 
model which demarcates the approximate locations of 
the subcritical and supercritical percolation phases. Al¬ 
though we haven’t attempted to determine the exact lo¬ 
cation of the line of criticality which separates these two 
phases, it is clear that our system lies within the super¬ 
critical percolation phase. 

Figure [9jr shows the spanning probability we obtain 
along a one-parameter path through our configuration 
space. The path, parameterized by e, travels along 
Pi = P 2 = | e for 0 < e < 1. The marginal proba¬ 
bility of obtaining a single bond in our lattice is pb = 
Pi +P 2 ~ 2piP2 = e(l — ^e) along our path. A percolation 
phase transition is seen to occur at pb = 0.484 ± 0.001. 
For comparison, in Figure we show a simulation of 
the standard square lattice bond percolation problem, 
wherein bonds appear independently of each other with 
probability ps■ Using identical methods, we identify a 
phase transition at ps = 0.500 ± 0.001, in agreement 
with the known exact value of ps = / ■ 

These results, along with the percolation results in 
Figure 4, conclusively demonstrate the supercritical be¬ 
havior of the random graph states obtained in our state 
reduction protocol, thus proving our ability to prepare 
cluster regions within our Union Jack state using only 
Pauli basis measurements. 


2 . 

get 


Non-Clifford Gates using our Interaction Gad- 


We first prove that our interaction gadget, asso¬ 
ciated with the measurement pattern shown in Fig- 

r(123) _ 


10 1 , implements the 


unitary gate U\ 


ure 

CCZ^VCZ (12) VCZ (23) , and we give the Clifford 
byproduct operators associated with certain unintended 
measurement outcomes. We then discuss how such gad¬ 
gets can be embedded into a surrounding cluster region, 
allowing them to act on arbitrary triples of qubits within 
that region. 

The core of our interaction gadget is the three-body 
operation given by multiplying two overlapping copies of 
CCZ and contracting one of the overlapping sites with an 


a) 



b) 


Qi / 


Q3 


r( 123). 


Qi 


Q2 


FIG. 10. a) Our interaction gadget, which allows us to apply 
the gate Ui to logical information. Blue triangles here rep¬ 
resent CCZ gates involved in forming the Union Jack state 
which play nontrivial roles in preparing Ui. We measure one 
control site in the Y-basis, six logical sites in the X-basis, and 
one logical site (along with many surrounding control sites) 
in the Z-basis, then use postselection to fix 13 of the control 
site measurement outcomes. The postselection is necessary 
to guarantee we can teleport information through the interac¬ 
tion gadget, the teleportation being carried out with the six 
X-basis measurements, b) The three-body operation which 
produces the diagonal unitary gate, Ui. Qubit 4 is initialized 
in a |+X) state, then contracted with a (+Y| outcome. 


ancilla state |+A”) and a Y-basis measurement outcome 
(±Y| = ^=((0| q= i(l|). Choosing (+Y| to be the ideal 
outcome, this yields the operation 


I/j 123) = (+Y| 


(4) 


(CCZ^CCZ^ 234 ^ |+X) (4) , (D3) 


which is diagonal in the computational basis (shown 
in Figure 10 3 ). Up to overall normalization and 


phase, 

on 1110) 
otherwise 


(1237 


Ui 

(123) 


ccz^Vcz il 2 ) Vcz {23) 


gives a phase factor of i when acting 
)r |011) (123) , and a phase factor of 1 

proving that its operation is identical to 
_ 

Because (— Y\ = ((+F|)*, 
the operation given by the outcome (— Y\ is (C/j 123 ^)*, 

(' 123 ) 

which is equal to U\ ' up to Clifford byproduct opera¬ 
tors czwczw. 

The three-body operation discussed above assumes 
that a (0| outcome has been obtained in the logical site 
Z measurement adjoining the control site Y measure¬ 
ment (yellow Z in Figure 10 1 ), and thus needs to be 
modified when a (1| outcome is obtained. In this lat ter 
case, the overlapping CCZ^^CCZ^ 23 ^ in Eq. (D3) is 
replaced by CCZ^CCZ^CZ^CZ^, and it can 

( 123) 

be shown that the resultant gate is again equal to U] 
up to Clifford byproduct operators S^ 3 \ where S is 
the phase gate S = diag(l,z). Finally, the case of unin¬ 
tended Y and Z outcomes in conjunction leads to Clifford 
byproduct operators (CZ^CZ^ S< 3 >)t. 

In summary, we have shown that a combined Y and 
Z measurement is capable of converting two non-Clifford 
CCZ gates into a three-body non-Clifford 17/ gate, with 
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FIG. 11. a) A pattern of control site outcomes which pos¬ 
sesses the “correct wiring” for our interaction gadget (X-basis 
measurements not shown). The wires percolate towards sep¬ 
arate points on the boundary without intersecting each other 
and without being acted on by stray CZ gates, b) An in¬ 
correct wiring pattern, which would require us to try again 
somewhere else in order to obtain a usable interaction gadget. 
Note that such regions can still be used as portions of cluster 
regions, without impacting the overall efficiency of our proto¬ 
col. In this case the control site marked Y would instead be 
measured in the computational basis. 


variation in measurement outcomes being accounted for 
by Clifford byproduct operators. Now how do we use this 
three-body unitary as a logical operation? One method 
for doing this is by measuring the control sites surround¬ 
ing our interaction gadget in the computational basis, 
and then attempting to use the random graph state we 


obtain to teleport qubits through the sites which Uj acts 
on. In the process of teleporting this information, Uj is 
successfully applied to the three qubits of interest. How¬ 
ever, we aren’t guaranteed to obtain a graph state with 
the “correct wiring”, i.e. one for which we can separately 
teleport each logical qubit to and from its respective site 
adjoining the interaction gadget, as in Figure |TT[ i. Be¬ 
cause of the possibility of obtaining graph states with 
incorrect wiring patterns, the successful embedding of an 
interaction gadget into a surrounding cluster region only 
occurs with some probability, which generically depends 
on the size of the surrounding cluster region. 

We can show that the probability of obtaining a cor¬ 
rect wiring pattern in the large system limit is finite and 
non-zero, by exploiting the same supercritical percola¬ 
tion properties which allowed us to prove the successful 
preparation cluster regions. This constant success prob¬ 
ability then guarantees that the stochastic nature of our 
interaction gadget embedding only contributes a constant 
multiplicative factor to the number of sites measured in 
our protocol. Consequently, our MQC protocol gives a 
proof of principle demonstration that we can efficiently 
perform quantum computation. Our proof involves first 
restricting ourselves to a region of finite size surrounding 
a particular interaction gadget, then using postselection 
(with finite success probability) to obtain a pattern of 
control qubit measurement outcomes which prepares a 
graph state with the correct wiring. For example, choos¬ 
ing a 6 x 6 grid of control qubits, we could postselect for 
the pattern shown in Figure [lTJi. 

When our region is of sufficient size, our postselected 
pattern can always be chosen so that distinct logical wires 
percolate without intersecting each other, and end at suf¬ 
ficiently separated points on the boundary of this region. 
When the separation between adjacent wire endpoints on 
the boundary of our finite region is much greater than 
the characteristic percolation length scale (the length 
scale associated with the exponential decay seen in Fig¬ 
ure 4), the conditional probability of continuing our post- 
selected pattern to a macroscopic graph state with the 
correct wiring factorizes into six uncorrelated probabili¬ 
ties. These probabilities, one for each wire, encode the 
chance of each wire percolating to a point infinitely far 
from its starting point on the finite region boundary. Be¬ 
cause of the supercritical nature of this percolation, each 
of these conditional success probabilities is finite, mean¬ 
ing that the total success probability for embedding an 
interaction gadget in a large cluster region is finite. Thus, 
our interaction gadgets can be embedded in cluster re¬ 
gions with a constant multiplicative overhead, letting us 
efficiently use them as logical gates which, together with 
the Clifford gates we get from our cluster regions, form 
a universal gate set. 
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